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Abstract 



We investigate what conditions allow us to extract the relative distribution 
of freeze-out space and time points in an arbitrary reference frame using 
tomography and source imaging. The source function may be extracted from 
the two-particle correlation function measured in heavy-ion collisions using 
imaging techniques. This imaged source function is related to the relative 
distribution of freeze-out space and time points through a generalization of 
the Radon transform found in tomography. Using tomography, the imaged 
source function may be converted into the relative freeze-out distribution in 
the frame of interest. We describe how the tomography may be performed in 
practice. 

PACS numbers: 25.75.Gz, 25.75.-q, 42.30.Wb 
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Typeset using REVT[t;X 



We intend to create the quark-gluon plasma in central nuclear collisions at the Relativistic 
Heavy-Ion Collider (RHIC). The quark-gluon plasma is predicted to have a long lifetime [jl],|^ 
which should lead to large time separations between the emission (or freeze-out) times of 
pairs of like particles. As two-pion correlations are sensitive to this time separation, they 
should be a useful tool for studying the plasma ||^,|^ . The signal of a long lifetime would be 
anomalously large correlation radii, particularly in the longitudinal and outward directions 
(in a analysis done in the Bertsch-Pratt parameterization |^). The lifetime information is 
encoded in these directions in a non-trivial manner which is currently accessed only using 
model parameterizations of two-pion correlation functions ||^. A long lifetime should also 
effect other like-pair correlations, such as two-proton or two-kaon correlations • However 
as in the case of protons, final state interactions often obscure this information. 

Recently it was shown that one can perform model-independent extractions of the en- 
tire source function S'(r) from two-particle correlations, not just its radii, using imaging 
techniques P-p!T|. Furthermore, one can do this even with relatively complicated final-state 
interactions (e.g. protons) and without making any a priori assumptions about the source 
geometry or lifetime, etc. First results from the application of imaging to the proton, pion, 
kaon, and IMF correlation data can be found in Refs. . These results are both intrigu- 

ing and nearly as hard to interpret as the original correlation functions: the reconstructed 
sources are found in the pair Center-of-Mass (CM) frame and the time dependence is folded 
into the sources in a non-trivial way. 

These two ambiguities are intimately tied together: the time direction is folded into the 
source function by a line integral in the direction of the boost from the system (or lab) frame 
to the CM frame. As we show in Eq. (Q), this line integral is a generalized Radon transform 
of the relative distribution of space-time emission points. If we construct the correlation 
function from pairs with a fixed rapidity and transverse momentum, then we know the 
boost from the system frame to the pair CM frame. Therefore, it should be possible to use 
tomography to reconstruct the relative distribution of space-time emission points from the 
imaged sources and we give an explicit inversion formula in Eq. ([T3|). There is one caveat: 
for the reconstruction to be unique, we must require that the probability for emitting a pair 
with a certain relative separation be independent of the total pair momentum. In what 
follows, we will refer to this condition as the momentum averaging approximation. The 
experimental signal for this condition to hold would be a sideward radius parameter that 
does not depend on the total pair momentum and this condition seems to be fulfilled at the 
CERN-SPS |T^. This analysis may be useful in the correlation function analysis of other 
systems, such as fragmenting hadronic strings or relativistic atomic collisions. 

The outline of this letter is as follows. First, we detail the imaging of the source function 
as a function of pair momentum and show that the imaged source function is a generalized 
Radon transform of the full two-particle source. Second, we study the reconstruction of 
the relative space-time distribution of emission points from the imaged source. Finally, we 
outline how to perform this reconstruction in a practical manner. 

We begin with a discussion of the correlation function and its relation to the source 
function. For concreteness, we consider like-pairs from central relativistic nucleus-nucleus 
reactions. We are eventually interested in quantities in the lab frame, namely the center-of- 
mass frame of the colliding nuclear system. The pairs we consider can be protons, pions, etc., 
as the formalism works for any like-pairs. The like-pair correlation function in an arbitrary 



2 



frame is the following ratio of invariant spectra: 



JTytruc / JTymixed 

C(P, q) = E,E,^l— /e,E,^^1—- . (1) 
apidp2 I apidp2 

Here, dN^^^^° / dpidp2 is the two-particle spectrum of pairs from the same event (averaged 
over many events) and dN'^'^^'^'^ / dpidp2 is the two-particle spectrum constructed from pairs 
in different events. In this expression, pi and p2 are the three-momentum of each of the 
particles and Ei and E2 are the on-shell energies. For simplicity, we work in side-out-long 
coordinates (a.k.a. Bertsch-Pratt [Q) coordinates where the longitudinal axis is along the 
beam line, the outward axis is along the component of P perpendicular to the longitudinal 
direction, and the sideward axis is perpendicular to the other two directions. In these 
coordinates, the total four-momentum of the pair is P = + p2 = {E, Pj^, Pq, 0) and the 
relative four-momentum is g = ^(pi — P2) = {qo, qi, lo, Qs)- 

In the Koonin-Pratt formalism, the correlation in Eq. (|1|) is related to the normalized 
single-particle sources through [p|,p|, [l0| , p!4| : 



C(P,g)= J d^r\<l>'^-^'^\r,q)\^ J d^R D{R + r /2, P/2 + q)D{R - r /2, P/2 - q) (2) 

Here $(~)''''^(r, q) is the pair relative wavefunction of the emitted pair and it includes all final- 
state interactions between the pair as well as (anti-)symmetrization due to statistics. Also 
D{r,p) is the normalized single-particle source and it denotes the probability of emitting an 
on-shell particle with four-momentum p at position f at time t. In terms of the emission 
rate, D{r,p) is: 

, Ed'^N /Ed^N 
d: p d T I d p 

where the emission rate is the number of particles frozen out per unit time, per unit volume, 
per unit volume of invariant momentum. Note that the normalized single-particle source 
transforms as a four-scalar. 

We may measure the correlation function and write Eq. in whatever frame we wish. 
In the pair center-of-mass (CM) frame Eq. (Q) simplifies considerably: 

Cp{q') = I rfV |$(-)-i(r',gOpy' d^R' DpiR' + r' /2,q')DpiR' - r' /2, -q'). (4) 

The subscript P refers to the momentum of the boost from some other frame (e.g. the lab 
frame) to the pair CM frame and the primed coordinates are pair CM coordinates. Note that 
in the pair CM frame P' = and E' = 2m. The reader should also note that each pair has a 
different P so it is boosted to a different pair CM frame. Furthermore, because of our choice 
of coordinates, this boost is confined to the longitudinal and outward directions. In order 
to remove the t' dependence of the wavefunction from Eq. (^), we have taken advantage of 
the facts that the particles are on-shell, giving qo = (3 ■ q, and that our wavefunction is a 
function of the Lorentz scalar q ■ r = {tj3 — f) ■ q = —r' ■ q'. 

We now make two crucial assumptions: the smoothness approximation and a momentum 
averaging approximation. For the smoothness assumption, we assume that the normalized 
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single-particle sources have only a weak dependence on q' . This is justified because the 
single-particle source varies weakly in r' as a function of q' while the wavefunction oscillates 
rapidly in r' for large q' ||T^ , p!5[| . Thus, at low q' the single-particle source is approximately 
independent of q' while at large q' the integral of the source function and the wavefunction 
average to zero. This approximation allows us to define the source function PJT0|1: 



Sp{?) = J dt' J d^R! Dp{R' + r'/2, q' ^ 0)Dp{R' - r'/2, q'^0). (5) 

Sp{r') is the probability of creating a pair a distance r' apart in the pair CM frame. In 
terms of the source function, Eq. becomes 

Cp{q')=[ c/V |<l>(-><=i(r',(p)p5^(r'). (6) 



Since the source function in Eq. (j^) is independent of the relative pair momentum, we 
can uniquely invert Cp{q') to obtain Sp{r') using the imaging techniques of Refs. P,p!0 



Therefore, we can dispense with the correlation function and work directly with the imaged 
source function. 

For reasons that will become clear momentarily, let us define the two-particle source as 
the probability of creating a pair of particles space-time distance of r = (t, r) apart in the 
lab frame. We may write the the two-particle source 5(r, P) in terms of the normalized 
single-particle sources as 

S{r,P) = I d'^RD{R + r/2,P/2)D{R-r/2,P/2) 

- ~ ^ ~ (7) 

d^R' Dp{R! + r72, q' ^ 0)Dp{R' - r'/2, q'^0). 



The second line in this equation is just the first line rewritten in the pair CM frame. The 
space-time displacement r' in the second line should be understood as a function of the boost 
velocity and space-time displacement r in the original frame. Here we have dropped the q 
dependence in both the lab frame (as in the lab frame P ^ q) and in the pair CM frame 
(in the smoothing approximation). Comparing the definition of the source function (P with 
the second line of Eq. (|^), we see 



Spir' 



J dt'S{r,P). (8) 



Thus, the imaged source is a line integral of the two-particle source along some boosted 
direction given by the pair total momentum, P. 

In the momentum averaging approximation, we replace S{r, P) in Eq. (^ with S{r). Here 
the S{r) has the interpretation as the relative distribution of space-time emission points, i.e. 
it gives the probability density for emitting any pair from space-time points separated by a 
displacement r in the lab frame. In terms of this distribution, we write Eq. (^) as 

Sp{r') = I dt'Sir). (9) 

Here all of the P dependence of the source function is attributed to the boost from the lab 
frame to the pair CM frame. This transform may be regarded as a generalized Radon trans- 



form as it involves boosts instead of rotations. In the conventional Radon transform |[T7|-[19 
one integrates a function /(x, y) along a direction rotated at an angle 9, e.g. 
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feix') = j dy'fix,y) (10) 

where y' = x sin 9 + y cos 9 and x' = x cos 9 — y sin 9. In Eq. (^, we are integrating along a 
time direction in a boosted frame. Thus, t' in Eq. (|^) plays the role of x in the conventional 
Radon transform and P (or equivalently the boost velocity) plays the role of the rotation 
angle 9. 

The momentum averaging approximation is justified for the following reasons. First, we 
could just define the relative distribution of emission points through Eq. (^. In this case, we 
will be folding any P dependence (such as from collective motion) into the time direction. 
Alternatively, we may think of S{r) as 5(r, P), but averaged over P. This is the viewpoint 
we adopt in this letter. Finally, in some systems or situations the two particle source only has 
a weak dependence on P and we may neglect this dependence. In a heavy-ion reaction, we 
can check this approximation by measuring the (e.g. pion) correlation function as a function 
of total pair momentum. Then we must examine the dependence of the correlation in the 
sideward and longitudinal/outward directions. If the dependence of the correlation in the 
sideward direction is relatively flat as a function of the total pair momentum, the two particle 
source does not have a P dependence. We will illustrate this below in a simple Gaussian 
model. We comment that the sideward pion radius parameter measured by the CERN- 



SPS experiment NA49 is relatively flat as a function of pair rapidity and transverse 
momentum so the momentum averaging approximation seems justifled. Strictly speaking, 
comparing the outward and sideward direction tells us that there is no P dependence of 
S{r,P) in the sideward/outward directions, but tells us nothing about the P dependence 
in the longitudinal direction. We postpone the detailed investigation of this approximation 



(especially in the presence of collective motion) for a future article ^0 . 

Let us now illustrate how the transform in Eq. (P) works on a model two-particle source 
with time evolution, but no dependence on the total pair momentum. Assume that we have 
spherical stationary two-particle source with radius R and time duration r in the lab frame: 

S«,f)=S„exp(-|5-^). (11) 

Since S{t, r) denotes the probability density for producing the pair with a space-time sepa- 
ration (t, r), the normalization constant is Sq = {Att'^tR'^)^^ . The boost from the lab to the 
pair CM frame is characterized by the boost velocity f3, which may be written in terms of 
the total pair momentum giving (3 = P/E = {(3l, (3o-,Q)- The coordinates transform from 
the pair CM frame to the system frame via f = r' + -^{P ■ '"O/^ + iP'^' ^ — + /3 ■ r') 

with 7 = 1/ y'l — P"^. Integrating over the CM time to flnd the imaged source function, we 
flnd: 




Sgir') = So^,i-^^,e^^ ( --4, K - ) ) ■ (12) 




This source function is also Gaussian and we immediately see how the radius parame- 
ters are modifled by the boost. First, since (3 has no component in the sideward di- 
rection, the Gaussian radius in the sideward direction is independent of (3 and is equal 
to the radius of the original two-particle source, Rs = R- We remind the reader that 
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a flat sideward radius parameter is a signal of the validity of the momentum averaging 
approximation. Second, both the outward and longitudinal radius parameters become 
Ro/L = R/ ^1 — Po/l(^'^ + t'^)/{E? + /J^r^). So for large emission durations, we flnd 
anomalously large radius parameters in the outward and longitudinal directions. Finally, we 



flnd a long/out cross term Rl^ = -R^PLpo^iR^ + r^)/{R^ + P^r^)- 

Now we can invert Eq. (§) and obtain the inverse in the Bertsch-Pratt coordinates: 

S{t,rL,ro,rs) = 

^ J dr', dr'o J dp, dPo g {P ■ {P -^{f-tp)))^S^^{r'^,r'o,r's). ^^^^ 

Here g{z) = dfi /x^e~*'^^ is a universal fllter function. Due to our choice of coordinates, 
rs = r'g. The derivation of the inverse in ([13|) is straightforward following the appendix of 



Ref. [|T6|, so we postpone the detailed derivation for a future article pO[| . 

Let us now comment on inverting Eq. (P) in practice. We might just plug the imaged 
source into Eq. (0) and integrate to flnd the relative distribution of space-time points. This 
method is called the flltered back-projection algorithm and it is a poor way to perform the 
inversion in practice because of the singular behavior of the fllter function g. A regularization 



of the fllter function is possible, but introduces numerical inaccuracy |]T8[. Furthermore, the 
filtered back-projection algorithm is computationally intensive owing to the four integrals 
in Eq. (M). 



A better approach is to use Algebraic Reconstructive Tomography (ART) [p!8i|T9| . In 
ART, one expands both the imaged source function and the relative distribution of space- 
time emission points in Eq. (P) in some function basis. For the sake of illustration we 



will use Basis splines ||2T| for this basis. In this basis, Sp{r') = J^ij SijBi{P)Bj{r') and 

S{r) = J2k^kBk{r) where and k are indices in the directions of P,r', and r respectively. 
Eq. (^ may then be recast as 

Y,S,,B,{P)Bj{r') = J dt' Y^SMr). (14) 

ij k 

Since the Basis splines are orthogonal, one can write this as a matrix equation: 

S^, = Y.^k J d^P d'r' B,{P)Bji?) J dt'Bkir) (15) 

k 

Only the t' integral needs to be done numerically as the other six integrals may be done an- 
alytically. Thus, ART reduces to inverting a matrix equation for the coefficients of the Basis 



spline expansion. We will explore the viability of the this approach in a future article pO |. 

In conclusion, we should be able to reconstruct the relative distribution of emission 
space and time points from the correlation functions measured in heavy-ion collisions. Our 
approach is model-independent and works for any like-pair correlations as it is implemented 
at the level of source functions. The momentum averaging approximation involved in the 
reconstruction seems to be justified in some experiments and its validity can be checked 
experimentally. Thus, this reconstruction provides us with a way to directly measure the 
freeze-out duration of heavy-ion collisions. 
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